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ABSTRACT 





The identification of the subsoil constitutive materials, as well as the detection of possible interfaces and anom- 
alies, are crucial for many site characterization applications. During investigation campaigns, complementary 
geophysical and geotechnical methods are usually used. These two sets of methods yield data with very different 
spatial scales and different levels of incompleteness, uncertainty and inaccuracy. In this work, a mathematical 
combination of geophysical and geotechnical information is proposed in order to produce a better subsoil char- 
acterization. It is shown that belief functions can be used for such a fusion process. A specific methodology is de- 
veloped in order to manage conflictual information and different levels of uncertainties and inaccuracies from 
different investigation methods. In order to test and validate this methodology, we focus on the use of two se- 
lected methods, Electrical Resistivity Tomography (ERT) and Cone Penetration Test. First, a synthetic model 
with artificial data is considered, taking advantage of the results obtained to conduct a comparative study (effect 
of parameters and noise level). Then, an experimental test bench is considered, in which a two-layered model is 
placed (plaster and saturated sands) and geophysical and geotechnical data are generated, using a mini-ERT de- 
vice and insertion depth values. This work also aims at providing a better graphical representation of a subsoil 
section with associated degrees of belief. The results highlight the ability of this fusion methodology to correctly 
characterize the considered materials as well as to specify the positions of the interfaces (both vertical and hor- 


izontal) and the associated levels of confidence. 


1. Introduction 


For subsoils characterization, investigation campaigns are set up, 
usually consisting of geophysical and geotechnical methods. These 
two families of methods are complementary and are used for various is- 
sues such as the characterization of slope stability (Caris and Van Asch, 
1991; Merritt et al., 2014; Chambers et al., 2013; Abidin et al., 2012) the 
characterization of potentially dangerous sites (James et al., 2014), the 
characterization of sites at construction (Coker, 2015) or the character- 
ization of river embankments (Perri et al., 2014). 

On the one hand, geophysical methods are non-intrusive and pro- 
vide physical information on large volumes of soils but with significant 
potential uncertainties. These uncertainties are due in particular to the 
integrative and indirect aspects of the methods as well as to the resolu- 
tion of the inverse problems. On the other hand, the geotechnical inves- 
tigation methods are intrusive and provide more punctual information 
but also more accurate. An important issue for the assessment of 
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subsoils is to be able to combine acquired geophysical and geotechnical 
data, while taking into account their respective uncertainties, inaccura- 
cies and spatial distributions (Royet et al., 2013). The complementarity 
of these two sets of methods is often underused since the uncertainty 
and inaccuracy associated with each method are rarely considered. Fur- 
thermore, the results are usually only graphically superimposed 
(Shaaban et al., 2013) instead of being mathematically merged. 

To characterize a section of subsoil and its potentially risky areas, it is 
essential to distinguish the different materials in place. The horizontal 
and vertical interfaces, as well as possible anomalies, have to be located. 
For levee embankment, as an example, it is in these locations that inter- 
nal erosion is likely to develop, which may lead to the complete rupture 
of the levee (Foster et al., 2000). Such a section characterization, with 
associated confidence indexes, could be included in failure hazard 
models. 

The use of belief functions (Shafer, 1976; Dempster, 1967) and dif- 
ferent information combination rules to combine geotechnical and geo- 
physical data is proposed. This makes it possible to take into account at 
the same time the uncertainties, inaccuracies and incompleteness of 
data related to each method. In the field of geosciences, belief functions 
have already been used and provide interesting results for slope 
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instability mapping (Binaghi et al., 1998; Althuwaynee et al., 2012), de- 
tection of precious metal (Tangestani and Moore, 2002), groundwater 
(Mogaji et al., 2015) or flood susceptibility mapping (Tehrany and 
Kumar, 2018). To our best knowledge, no work has been proposed, con- 
sidering the combination of two sources of information with different 
Spatial distribution (spatialized and punctual) and for an investigation 
campaign in the vertical section. 

Here, an innovative method of information fusion to combine elec- 
trical resistivity tomography results and cone penetrometer test data 
is proposed. First, work on data obtained from synthetic models is 
displayed. The obtained results allow to conduct a comparative study, 
evaluating the effect of different parameters (like the data noise level) 
on the fusion result. The fusion methodology is then tested from data ac- 
quired on a test bench. In this work, the potential of such a methodology 
is shown by using insertion depth data, acquired by a laboratory pene- 
tration cone, and electrical resistivity data acquired by a mini Electrical 
Resistivity Tomography (ERT) device. The depth of penetration data 
corresponds to geotechnical information while the electrical resistivity 
data correspond to geophysical information. The main concern is to 
highlight the ability of this information fusion algorithm to characterize 
the interfaces between materials and to discriminate three different 
types of materials with variation in thickness of one of them, and to 
present the variation of the results according to the number and posi- 
tion of the simulated boreholes. 

The main contributions of this work are as follows. First, this new 
methodology makes it possible to take into account the uncertainties, 
inaccuracies and incompleteness associated with the different methods 
of investigation used, proposing a modeling of the Basic Belief Assign- 
ments (BBAs) specifically adapted to the problematic. Then, the pro- 
posed graphic representation is innovative since it allows both to 
present the different geological sets that would be present in the subsoil 
and their layout, while presenting the confidence associated with these 
results. This methodology is particularly suitable for the characteriza- 
tion of interfaces and anomalous zones, which may correspond to 
areas where the risk of instability is potentially the greatest. This work 
also allows the implementation of a small physical model to validate 
the fusion approach with real data. 

This article is organized as follows. In Section 2 a presentation of the 
approach of fusion used in the methodology is given, which introduces 
the use of the evidence theory and the combination methods used here. 
In Section 3, a synthetic study will then present the fusion approach 
from artificial data. It will also present the comparative results associ- 
ated to two parametric studies. Then, in Section 4, a presentation of 
the investigation methods used in the introduced experiment (labora- 
tory penetration cone and mini ERT device) is given. Finally, the test 
bench fusion results are presented in Section 5 and discussed in 
Section 6, in order to understand the interests, limitations and perspec- 
tives of such a methodology. 


2. Fusion methodology 
2.1. Belief functions and combination rules 


The belief functions have been introduced by Shafer (Shafer, 1976) 
in 1976 in the development of his mathematical theory of evidence in- 
spired by previous works of Dempster (Dempster, 1967). Shafer's the- 
ory is also referred as Dempster-Shafer theory (DST) in the literature. 
This theory (proposes a method to) calculate(s) the belief and the plau- 
sibility of an event (here a soil material class) from distinct source of ev- 
idence (measured data). The practical advantage of using such a theory 
lies in its ability to manage information from different sources, associ- 
ated with variable uncertainties and inaccuracies. In this work, only 
two sources of information will be considered: geotechnical and geo- 
physical. Another advantage of this theory is its ability to assess the de- 
gree of conflict between sources (ex: contradictory information 
between data obtained from large scale geophysical campaign and 


from punctual geotechnical investigation). Uncertainties correspond to 
degrees of confidence that are given to a value, whereas inaccuracies 
correspond to intervals of values that can be directly associated with 
measurement errors related to the method. For example, the uncer- 
tainty of measuring a geotechnical parameter identical to the one mea- 
sured in a borehole increases with the distance to that point. The 
inaccuracy can for its part, be associated with the error bar of the result. 
The belief functions allow to take into account the ignorance and incom- 
pleteness of the information. It is indeed possible to grant credit on all 
the possible results in order to quantify the ignorance. For the reader 
eager to learn more, the theory is detailed in (Martin et al., 2008). 

A Bayesian approach as part of a subjective probability approach 
(Cooke, 1991) could have been considered for geophysical and geotech- 
nical data combination. However, the main limitation of such an ap- 
proach is that probabilities essentially represent uncertainty and only 
very poorly the level of inaccuracy. Moreover, in the probabilistic 
modeling stage, the different decisions (events) are only represented 
on singletons (i.e. single events ) and are necessarily considered exhaus- 
tive and exclusive. The exclusivity is implied by the assumption of the 
additivity of probabilities. However, this hypothesis may be too strong 
and limit the representation of the knowledge. Furthermore, with a 
Bayesian approach, it is difficult to model the lack of knowledge or the 
knowledge that is not expressed in probability distributions. 

In order to define and to use the belief functions, it is necessary (1) to 
set a frame of discernment, (ii) to assign belief mass values to the events 
of this framework (Basic Belief Assignments - BBAs), (iii) to choose a fu- 
sion rule for combining information; and (iv) to represent the combined 
information. 

The Frame of Discernment O (FoD) is made of all the possible events 
about the problem under concern, the elements of the FoD are exclusive 
and exhaustive, so that for n events: 


© = {6),62,...,On} (1) 


In the considered problematic, the possible events of the FoD corre- 
spond to intervals of values of geophysical and geotechnical parameters 
that can be associated with classes of geological materials (for example, 
0, = clays, 02 = sands...). The space of belief mass functions, the set of 
all subsets of O, written 2°, is fixed by all the disjunctions and by the 
possible conflict between the sources of information (written @) such 
that: 


2° = {Ø, {01}, {02}, {01U82}, {03}, {01083}, 2) 
{02U63 }, {8;U02U03}, ..., {0; U82U03U...U8n }} 

As in the probability theory, the belief mass function mj is defined, 
for a source of evidence S; (for j = 1, 2), attributed to A (defined on 
2°) in [0, 1] such that the more m(A) tends to 1 and the more the confi- 
dence in A is important: 


>> mj(A) = 1 (3) 


Ac2® 


The difference with the probability theory lies in the fact that A can 
represent the union of several events (for example, either 0; OR 62). It 
is therefore possible to model uncertainty and lack of knowledge. For in- 
stance, when no information is available about the achievement of an 
event member of O, one can set m;(@) = 1, avoiding the uniform distri- 
bution that would have been considered in a probabilistic scheme. Com- 
bination rules, as part of the belief functions theory, can thus take 
different levels of uncertainty and imprecisions into account according 
to the source of information. If only defined on singletons, the belief 
mass function is similar to a probability distribution. 

Smets fusion approach developed in his Transferable Belief Model 
(TBM) (Smets, 1990) (i.e. conjunctive fusion) allows the attribution of 
a mass of belief to the conflict, outside the FoD, so that (open-world 
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assumption): 
m12 ()>0 (4) 


where ™m2(-) denotes the combined BBA resulting from the combina- 
tion of information of sources 1 and 2. The belief mass resulting from 
the fusion of information from source 1 and 2 is written: 


mı2 (A) = mı(X)m2 (Y) (5) 
X YCOIXNY=A 


And the level of conflict between the two considered sources of in- 
formation can therefore be quantified by: 


m12 (Ø) = mı(X)m2(Y) (6) 
X YCOXNY-Ø 


With m,(X) and m2(Y) the belief masses respectively attributed to 
events X and Y by sources 1 and 2. 

According to Shafer's approach and unlike Smets' rule, Dempster- 
Shafer's rule (DS) does not allow the attribution of a mass of belief to 
the conflict (closed-world assumption): 


m3 (Ø) = 0 (7) 


The conflict is there reallocated through a classical normalization 
factor. The mass of belief in A, m73(A), resulting from the fusion of infor- 
mation from sources 1 and 2 is written: 


1 


m3 (A) = TaD 


mı(X)m2(Y) (8) 
X,YCOJXNY=A 





The disadvantage of this method is that the conflict between the 
sources is no longer represented and it is possible to obtain counterintu- 
itive results if the conflict is important because of this normalization. 
Even more problematic, even if the distinct sources are both informative 
whatever the level of conflict is, Dempster-Shafer's fusion process can 
even not take into account the second source of information (Dezert 
et al., 2012). 

On the other hand, the PCR6 (Proportional Conflict Redistrib- 
ution No. 6) rule of combination (Smarandache and Dezert, 2009; 
Smarandache and Dezert, 2005), considering two sources of informa- 
tion, only transfers the conflicting mass to the events that are actually 
implied in the conflict and in proportion with their individual masses 
in order to preserve the specificity of the information. In order to 
apply the PCR6 rule: i) the combination rule described in Eq. (5) must 
be applied, ii) the total or partial conflicting masses have to be 


0, 9, 05 


Number of cells 





a, a> a «ay as aç 






calculated and iii) the total of partial conflicting masses have to be 
redistributed proportionally on non-empty sets. So, for m{'$*°(0) = 0 
and Y A E 2°\{@}: 





E mı (A) m(Y) | _m(A)*m(Y) 
m AmA > mAr mA) O 
ANY = Ø 


2.2. Construction of BBAs from geophysical and geotechnical data 


Belief masses have to be assigned to each considered event of the 
FoD, for both sources of information. The combination of the belief 
masses can only be initiated after this stage. In the following, the geo- 
physical source of information will be identified as source 1 and the geo- 
technical source of information as source 2. A 2D model assumption will 
be made, corresponding to the x and z spatial axes, since vertical sec- 
tions of subsoil are considered. 


2.2.1. Geophysical data 

The discretization of the considered subsoil section, as well as the 
depth of investigation and the resolution, depend on the acquisition 
method used (Kearey et al., 2013). It is the user who sets, using the in- 
version tool used, the shape and dimensions of the discretization grid 
used. It is about starting from this discretization and being able to asso- 
ciate for each cell, masses of beliefs for each event of the FoD. 

The constitutive classes of the FoD are also fixed at the end of the in- 
version process by the geophysicist, with the help of a representation of 
the distribution of the set of inverted geophysical values, in the form of 
modal classes (Fig. 1a). The representation in this form makes it possible 
to highlight the centers, minima and maxima of the events considered 
in order to be able to fix the bounds of the intervals associated with 
the events of the FoD. The number of cells of the subsoil section are rep- 
resented according to the geophysical parameter values. The infima and 
suprema must be fixed so that the intervals are of the same width in 
order to avoid the appearance of a bias when calculating Wasserstein 
distances (detailed under). To associate the belief masses with the FoD 
events, the intervals of inverted values of the physical parameter (in 
red, Fig. 1b) are considered. For some geophysical methods, these inter- 
vals can correspond to the value obtained at the end of the inversion 
with its associated inaccuracy. 

It is then necessary to associate belief mass values m,(.) correspond- 
ing to each element of 2°, for each cell of the inverted section. The 
masses are obtained from the calculation of Wasserstein distances 
(Tran and Duckstein, 2002), considering two geophysical intervals 
A = [aj, a2] and B = [bj, b2] with A and B belonging to R, A being the 





b) 
Infima b, b, Maxima 


Geophysical 
i parameter values 


Fig. 1. a) Model classes’ distribution of the geophysical parameter values from the considered subsoil section, allowing the selection of the geophysical classes in b). The red interval 
corresponds to an interval of inverted values, from one cell of a 2D section of subsoil, used for Wasserstein distances’ calculation. (For interpretation of the references to colour in this 


figure legend, the reader is referred to the web version of this article.) 
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interval corresponding to an event of the FoD and B being an interval of 
inverted values (Fig. 1b), Eq. (10): 


min =[52)- (5%) 05 [OGY + OE) 





This calculation estimates the distance between two intervals ac- 
cording to their size and the distance between them. The Wasserstein 
distances are calculated (using a logarithmic scale if the geophysical pa- 
rameter requires it) between the inverted values with estimated inac- 
curacies, and the intervals associated with each event, chosen by the 
geophysicist. Each cell is finally associated with a standardized BBA re- 
specting Eq. (3). This way, the more the distance of a geophysical inter- 
val resulting from the inversion is “close” to one event of the FoD, the 
more the mass of belief associated is important, and reciprocally. 


2.2.2. Geotechnical data 

For the geotechnical part, the information proposed during an inves- 
tigation campaign is spatially punctual (in the x-z plane) and often 
contained in vertical soundings made from the surface. It is about asso- 
ciating masses of belief with the different events of the FoD for each cell 
of the considered vertical soundings. For this, the values proposed at 
each depth are considered with the associated inaccuracy, correspond- 
ing to the measurement error that could be attributed to the measuring 
device (Fig. 2a). Thus, as for the geophysical part, intervals of values are 
obtained. 

The geotechnical mesh consisting of as many cells in depth as the 
number of geotechnical values (Fig. 2b) is generated. A mass of belief 
m>(.) = 1 is assigned, in the drilling points, to the events corresponding 
to the measured geotechnical parameter. A value of 1 is set since we are 
very confident in the information inside the boreholes unlike the 
Spatialized geophysical information. A new mesh is then constructed 
(Fig. 2c), according to the size and depth of the boreholes. In order to 
characterize the entire section of the model, as does the geophysical 
method, and to associate mass values to each new cell (BBA), an expo- 
nential lateral decay of the belief mass is imposed, from the drilling 
point to the nearest borehole so that the decay rate is a function of the 
values proposed by the nearby borehole. So that, for a specific depth, 


Eq. (11): 


BBA(x) = e~*°* BBA(0) (11) 


With x being the distance from the considered cell to the reference 
borehole (x = 0 in the borehole), k a decay factor fixed by the user to ad- 
just the lateral decay rate, BBA(x) the belief mass values assigned to each 
event of the FoD for a position x, with BBA(0) = 1. C, corresponds to the 
coefficient of variation expressed in Eq. (12), such as used in (Phoon and 


SD 1 SD 2 SD 1 
X b) 


a) 






Geotechnical parameter values Z(m) 


Z(m) 





Kulhawy, 1999): 






— ~ (Q— Qi)? 
E, _ mesh (12) 
Q 





where Q is the geotechnical value of the reference cell in the considered 
borehole and Q; the geotechnical value in the nearby borehole. For Fig. 
2b, Nmesh = 3 has been considered. If Nmesn = 5 or 7, the computation 
of the C, will take into account 5 or 7 cells in the nearby borehole. In- 
deed, for two consecutive boreholes with similar values, at similar 
depth, the decay of the confidence is slower than for two consecutive 
boreholes presenting radically different values. This decay of belief 
mass is carried out to the left and to the right, from each drilling. 

If, between two boreholes, the mass of belief associated with a hy- 
pothesis A is <1 (mp (A) < 1), then the remainder of mass to be allocated 
to satisfy Eq. (3), is reported on the proposition “any type of material” 
represented by the union of all events, such as Eq. (13): 


məz (01 U02U03U04) = 1—mz(A) (13) 


Considering n boreholes, from 1 to n from left to right: borehole 1 
cannot be compared to any borehole to its left neither can borehole n 
be compared to any borehole to its right. Indeed, for a given depth, an 
equal C, is considered for left and right directions, for boreholes located 
at the beginning and at the end of the section. 


2.3. Dimensioning of the mesh prior to the fusion 


Each source of information imposes its own mesh but in order to 
combine the belief masses from the geophysical information source 
(source 1) and the geotechnical source (source 2), it is necessary to 
have a common mesh containing, for each cell, the geophysical and geo- 
technical BBAs. In order to not alter the quality of the information, no in- 
terpolation is carried out. It is decided to superimpose the geophysical 
discretization grid resulting from the 2D inversion to the geotechnical 
division, depending on the number and the borehole positions. Thus, 
an irregular mesh is obtained but without any approximation (Fig. 3). 


3. Synthetic study 


Below, a synthetic study based on artificial data is proposed in order 
to test this new proposed methodology. It is the opportunity to show 
the impact of different levels of noise on the geophysical information 
as well as the influence of the lateral decay factor k (Eq. (11 )) on the re- 
sults of the fusion in order to be able to choose a value for the use of such 
a methodology from real data. 


Sese e-eeeEe SEEE 
= =z 





























Fig. 2. Construction of a geotechnical discretization mesh from two vertical boreholes acquisition (SD1 and SD2). a) Representation of the geotechnical values for SD1 and SD2 according to 
the depth. b) The boreholes are divided in cells associated with belief mass equal to 1 for the considered event. c) Construction of a new mesh according to the size and depth of the 


boreholes. 
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Fig. 3. Example of a geophysical mesh (in black) and a geotechnical mesh (in red) 
superimposed to propose a new irregular mesh to carry out the combination 
calculations and present the fusion results. (For interpretation of the references to 
colour in this figure legend, the reader is referred to the web version of this article.) 


3.1. Considered methods 


For this study, the electrical resistivity tomography (ERT) method 
stands for the geophysical information source and the Cone Penetrom- 
eter Test (CPT) method for the geotechnical information source. 

The basic principle of DC-resistivity methods consist in injecting an 
electric current of known intensity [A] by means of two “current” elec- 
trodes and measuring a voltage [V] between two “potential” electrodes. 
Depending on the electrode layout, the topography, the properties of 
the materials and their distribution, apparent resistivity values can be 
computed. The depth of investigation depends on the spacing of the 
electrodes, the configuration of the electrodes and the nature of the 
soil (Loke, 2011). By generalizing this principle, a two dimensional 
(2D) ERT consists in aligning a series of electrodes and acquiring a 
large number of measurements based on four electrodes configuration. 
The apparent resistivity data acquired are then inverted using an inver- 
sion code or software to reconstruct a complete 2D-section of electrical 
resistivity [Q.m]. Here the Res2Dinv software (ver 3.71.118) (Loke, 
2011) has been used. 

In order to obtain an artificial resistivity section of subsoil, a two 
steps procedure is followed. First, resistivity data are simulated using 
the Res2Dmod software (Loke, 2002), on the section that we want to 
consider. Second, apparent electrical resistivity values are inverted 
with Res2Dinv, considering a L1 norm (Loke et al., 2003) and an ex- 
tended model discretization, to obtain the synthetic inverted section 
of electrical resistivity. 

The CPT method consists in pushing rods into the soil, at a con- 
stant speed, with a conical tip at the end (ISO 22476-1, 2012, n.d.). 
This test is often used for the determination of the soils mechanical 
resistance properties. The two measured parameters are the tip re- 
sistance q: [MPa] and sleeve friction fs [MPa]. Although the method 
uses two parameters, only qe will be considered as the study 
parameter. 


xX Om 


5m 


17m 
Om 24m 


3.2. FoD and considered model 


For this synthetic study, a two-layer model is considered, composed 
of materials that can be likened to silts for the upper layer and clays for 
the underlying layer. The FoD therefore contains three material class hy- 
potheses, such as: 


O = {01,02,03} (14) 


With 6, the event corresponding to the clayey material, 62 to the silty 
material and 63 to unknown materials. The latter is associated with 
the union of the geophysical and geotechnical value ranges that do 
not correspond to those associated with 6; and 62. This event 63 allows 
us in a certain way, to quantify the lack of knowledge of the environ- 
ment since it does not include the two first sets. The construction of 
the BBAs then consists in associating the data of the two considered 
sources to the events of the FoD. Fig. 4 shows the two-layer model 
based on events from the FoD, used for this synthetic study. 


3.3. Construction of BBAs from geophysical and geotechnical data 


3.3.1. Geophysical data 

The electrical acquisition is simulated with a Wenner acquisition 
mode and with 96 electrodes interspaced from one meter. An electrical 
resistivity of 100 0.m is considered for the upper material and a resistiv- 
ity of 30 O.m for the underlying one (Palacky, 1987) (Fig. 4). Electrical 
acquisitions are simulated with different noise levels (5, 10 and 15%). 
The results of this inversion (with 10% noise, Fig. 5a) allow to highlight 
the presence of two layers but the interface between these two layers is 
not perfectly identified. A variation of thickness in the center of the 
model is visible. The interface is not straightforward and anomalies 
are present on the surface even though they are not part of the initial 
model. 

From these inversion results, it is possible to define the ranges of 
electrical resistivities that will be associated with the different events 
considered for the fusion process. A distribution in modal classes is 
used to visualize the number of cells, in the discretized section of the 
2D inversion, associated with specific range of resistivities (Fig. 5b). 
This distribution allows to highlight the two large material classes of 
the model. Thanks to it, the bounds of the considered events can thus 
be defined (in 0.m), so that the intervals have the same length (in log- 
arithmic scale): 


6, = [25; 45] 
62 = [83; 149.4] (15) 
6; = [13.89; 25[U]45; 83[U]149.4; 268.92] 


As explained in II.2., it is possible to associate belief masses with each 
cell of the mesh thanks to the values resulting from the inversion. 
As part of the construction of geophysical BBAs, the values presented 
Fig. 6 are obtained. This figure highlights the association of the values 
of Fig. 5a with the events of the FoD, Eq. (15). The presence of a top 
layer (02) and a base layer (01) can be detected (Fig. 6a). It appears 
that there is a variation in the thickness of the layers in the center of 
the model, but the interface is not well characterized. Moreover, the in- 
termediate values of electrical resistivity resulting from the inversion 
(Fig. 5a) between 62 and 04 layers induce the representation of a third 
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Fig. 4. Representation of the events of the FoD in the imposed model of subsoil of the synthetic study. 
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Fig. 5. a) Subsoil section displaying inverted electrical resistivity values from 10% noise data acquisition and b) model classes’ distribution of the cells presented in a), according to the 


electrical resistivity values (Q.m). 


material (03) which has no reality in the model that has been fixed. The 
belief masses are maximum when the resistivity values correspond to 
the center of the resistivity classes set for each event (Eq. (15)). 


3.3.2. Geotechnical data 

Concerning the source of geotechnical information, the simulation of 
four vertical CPT soundings inter spaced from 19 m is proposed (Fig. 7). 
20 cm wide and up to 15 m deep boreholes are considered, and a value 
of qc is recorded every 50 cm from the surface. An inaccuracy of 
107? MPa on the measurements is considered. For a fixed normalized 
friction ratio of 3%, a value of q: of 20 MPa is considered for the upper 
silty material and a value of 0.2 MPa for the underlying clay material, 
as proposed in the Robertson diagram (Robertson et al., 1986). 

In order not to have uniform values of q< for the materials and to try to 
represent the noisy reality of an acquisition in the field, values are drawn 
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following a normal distribution defined for each event. Mean q, values of 
0.2 and 20 MPa are respectively used to define the normal distributions 
of the material classes. Standard deviation values equal to 10% of the 
mean values are associated, echoing the 10% noise used for the geophys- 
ical data. Keeping the minimum and maximum values, these random 
draws, make it possible to define the limits, in MPa, of the intervals asso- 
ciated with the elements (i.e. material classes) of the FoD: 


6, = (0.14; 0.27] 
6. = [13.5; 23.5] 
63 = [0.1;0.14[U]0.27; 13.5[U]23.5; 100] 


(16) 


The minimum and maximum values are fixed at 0.1 and 100 respec- 
tively because they are the minimum and maximum values in Robe- 
rtson's diagram (Robertson et al., 1986). 
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Fig. 6. a) Representation of the event having the highest belief mass according to the BBA construction from geophysical data and b) the associated belief mass values, considering a 10% 


noise. The black line represent the position of the interface. 
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Fig. 7. 2D section of subsoil displaying true ER distribution with boreholes positions in black and associated tip resistance vertical profiles in white. 


There are two types of sounding results according to their position 
(Fig. 8). Once the values associated with the meshes of the sounding 
are obtained, it is possible to associate masses of belief to the whole sec- 
tion by extending the geotechnical information, as explained in II.2. In 
the framework of the construction of geotechnical BBAs and for k = 
0.1 (Eq. (11)), the obtained values are proposed in Fig. 9. This figure 
highlights the fact that the confidence is maximum in the soundings. 
This method allows us to characterize the material 0» on the first 5 m 
of the model and the material 6, from 10 to 15 m deep. 

The greater thickness of the material 62 in the center of the model is 
also well characterized. On the other hand, a great doubt appears in yel- 
low (Fig. 9a) in certain areas, not allowing the determination of a spe- 
cific material (04 U 62 U 03). For the model base area, this can be 
explained by the fact that the soundings stop at 15 m depth. Regarding 
the areas between 5 and 10 m to the right and left of the first and last 
sounding, this is related to the fact that for these two soundings, the 
closest soundings propose different values at the same depths, the con- 
fidence attributed to the presence of 6,therefore decreases very quickly 
laterally. This decay is also high towards the edges of the model because 
no other sounding is present at the ends to constrain the information. 


3.4. Effect of lateral decay factor and noise level on the fusion results 


We examine in this part the results of the fusion of belief masses 
established for the proposed synthetic model by varying the noise 
level of the geophysical information, as well as the value of the lateral 
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decay factor k (Eq. (11)) influencing the lateral decay rate of geotechni- 
cal information. 

Fig. 10 shows the fusion results with different values of k (1077, 
5.1072, 107 ', 5.107! and 1) for a simulated noise of 10% on the acquired 
geophysical information. Noise was set at 10% since the electrical resis- 
tivity classes of the FoD were defined from the modal classes of the 
inverted 10% noise image, Fig. 5b. For each value of k, Fig. 10a and b rep- 
resent the results obtained by Smets fusion whereas Fig. 10c and d rep- 
resent the results obtained by PCR6 fusion. While Fig. 10a and c show 
the material classes having the greatest mass of belief at the end of 
the fusion process, Fig. 10b and d correspond to the values of these re- 
spective belief masses, between 0 and 1. These figures, display the 
events (materials) potentially present within the section, as well as 
their attached level of confidence. 

The higher the value of k is, the higher the rate of confidence in geo- 
technical information is. This can be seen, for example, from the last 
borehole to the right end of the section (Fig. 10b) or between the 2nd 
and 3rd boreholes (Fig. 10d). The increase of k implies that for two 
soundings offering similar values at the same depth, the confidence as- 
sociated with the corresponding type of material will tend to decrease. 
On the other hand, for two soundings proposing different values at 
the same depth, the increase of k will hardly have any impact on the be- 
lief masses associated with the selected events (e.g. between 5 and 10 m 
of depth between the boreholes 1 and 2, Fig. 10d). 

With regard to the material classes identified after the fusion pro- 
cess, the more k increases, the more the quantity of conflict de- 
creases (in red, Fig. 10a). This is explained by the fact that when 
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Fig. 8. Examples of the two types of simulated soundings with tip resistance values according to the investigation depth. a) corresponds to borehole n°1 and 4 on Fig. 7 while 


b) corresponds to borehole n°2 and 3 on Fig. 7. 
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Fig. 9. a) Representation of the events having the highest belief mass according to the BBA construction from geotechnical data and b) the associated belief mass values. The borehole 
positions are in dashed lines while the black line represent the position of the interface. 
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Fig. 10. Representation of the events having the highest belief mass in a) and c) and their associated mass values in b) and d), considering a 10% noise. In i) k = 0.01, ii) k = 0.05, iii) k = 0.1, 
iv) k = 0.5 and v) k = 1. Figures on the left side are results for Smets fusion while figures on the right side are results for PCR6 fusion. The sounding positions are in dashed lines while the 
black line represent the position of the interface. 
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there is little trust in the geotechnical data, there is little conflict 
with the geophysical data. In the meantime, an increase in the pro- 
portion of 63 is observed (Fig. 10c) close to the interface. This obser- 
vation is explained by a larger mass attributed to the union of events 
and by geophysical data which propose intermediate values at the 
interface level. 

In the following of this work, an intermediate value of k will be 
retained, equal to 0.1. With such parameter value, a good confidence 
in information repeating between two successive soundings is ob- 
tained, but it also leaves room for doubt by having enough unknown 
material (03) at the interfaces. The obtained fusion results with dif- 
ferent noise levels added to the geophysical information (5, 10 and 
15%) are shown in Figs. 10iii and 11 with k = 0.1. 

The greater the amount of noise is, the less clear the interfaces pro- 
posed by the inversion are (Figs. 5a, 11.1.a, ii.a). A greater number of 
anomalies are also present when the noise level increases. The noise 
level finally impacts the level of inaccuracy associated with the geo- 
physical data used in the fusion process. Larger data inaccuracies induce 
wider value ranges considered for calculating Wasserstein distances, 
which in turn can bring to consider belief masses on more events of 
the FoD. 
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Since the classes associated with FoD elements were fixed from the 
values with 10% noise (Section 3.3 and Fig. 4b), it is “reasonable” to 
have a higher confidence (higher belief masses) on these results than 
on the results with 5% and 15% noise (Figs. 11c, e, 10iii.b, iii.d). The fu- 
sion process allows to override the noise effects, whether the noise 
level is 5 or 15%. This can be imputed to the computation of Wasserstein 
distances, taking into account the data inaccuracies and considering all 
geophysical classes. 


4. Setting up a test bench for real geotechnical and geophysical 
acquisitions 


4.1. Materials 


In order to be able to assess the validity of the developed fusion 
methodology, two methods of data acquisition were retained: (i) a 
mini-ERT device acting as the geophysical source of information and 
(ii) a laboratory penetration cone acting as the geotechnical source of 
information. Before setting up the test bench, it was necessary to select 
the materials that could be put in place in a tank in order to carry out the 
study. This selection implies that the materials used meet several 





Unit electrode spacing 1.00 m. 





11.2 % 


Depth Iteration 3 RMS error = 
0.0 16.0 32.0 


0.0] 
1.65 
3.05 


4.74 
6.79 
9.27 


10.7 
12.3 


14.0 
15.9 


Inverse Model Resistivity Section 


E a a a c ee l o ‘l a ‘U ‘M ‘M ‘M 
25.0 31.4 39.3 49.3 61.9 77.6 97.4 122 


Resistivity in ohm.m 











Unit electrode spacing 1.00 m. 


EE 
ee. 
| Je, 
| - 

1 













Fig. 11. a) Subsoil section displaying inverted electrical resistivity values from i) 5% noise and ii) 15% noise data acquisition. Representation of the events having the highest belief mass in 
b) and d) and their associated mass values in c) and e). Figures on the left side are results for Smets fusion while figures on the right side are results for PCR6 fusion. The sounding positions 


are in dashed lines while the black line represent the position of the interface. 
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conditions in order to validate the methodology: they must have 
(1) distinct electrical resistivity ranges, (ii) distinct penetration depths 
and (ili) a certain homogeneity in the space to limit uncontrolled anom- 
alous values. 


4.1.1. Mini ERT device 

Expressly for the purposes of this study, a mini ERT device (Fig. 12) 
has been set up. This device consists of 48 electrodes of 6 mm length, po- 
sitioned at regular intervals of one centimeter. It can be moved along the 
test bench to make multiple acquisitions and to cover a longer section. 


4.1.2. Laboratory penetration cone 

The laboratory penetration cone method is described in the French 
standard (NF P 94-052-1, 1995). It consists in measuring a penetration 
depth of a cone, in millimeters, subjected to its own weight (Fig. 13). 
The materials are tested individually, repeatedly, to determine an aver- 
age value and a standard deviation of penetration depth for each mate- 
rial. These values can be used later in the study to simulate different 
drilling positions within the test bench. This method can be likened to 
the CPT method which is one of the most popular in situ geotechnical 
tests. 


4.1.3. Test bench and used materials 

For the validation of the methodology, we wanted to build a test 
bench that could be easily set up and controlled, with two or three 
layers and variation of the interface positions. Fast-hardening natural 
fine-grained plaster as well as Hostun fine sand (Flavigny et al., 1990) 
are the retained constituents. These two materials meet the three condi- 
tions listed above. They were placed in a transparent PVC tank of 100 
x 30 x 17 cm? as shown in Fig. 14 with an underlying layer of 5 cm of 
plaster (setting time = 69 h) overlaid by a layer of 2.5 cm of water sat- 
urated sand. 

A formwork was made during the placement of the plaster so that a 
20 cm long anomaly could be inserted in. Saturated sand of 7.5 cm thick- 
ness is present instead of plaster. The contact between the materials and 
the bottom of the tank is at the origin of an interface that will be inter- 
esting to detect with the help of the methodology. 16 kg of plaster were 
mixed with 8 kg of water to obtain the material finally put in place. The 
electrical resistivity of the plaster was measured before and after the 
placement of the saturated sand to verify that the presence of water 
had a negligible impact on the electrical properties of the plaster. 

For the Hostun sand, 15.82 kg were pluviated in 5.8 kg of water, 
above the plaster to reach saturation. Trials had been carried out in ad- 
vance to determine the proportions of water and sand required to 
achieve such a state as well as to validate the repeatability of such instal- 
lation by pluviation. The values of electrical resistivities and penetration 
depths are displayed in Table 1. 





Fig. 12. Mini ERT device with 48 electrodes spaced 1 cm, adjustable height, used for 
electrical acquisitions in the test bench. 





Fig. 13. Laboratory penetration cone. 


4.2. FoD and BBA modeling 


4.2.1. FoD and target model 
A FoD consisting of four elements (material classes) is considered so 
that: 


O = {01, 02, 03,04} (17) 


With 6, the element corresponding to the plaster material; 02 corre- 
sponding to saturated sand; 63 corresponding to the hard and electri- 
cally insulating bottom of tank simulating a substrate and 6, 
corresponding to unknown materials, being the union of the ranges of 
values not corresponding to those associated with the 3 previously de- 
scribed materials. Fig. 15 presents the target model in the form of events 
constituting the FoD, following the disposition of the materials within 
the test bench. Although the tank used is 1 m long, the ERT acquisition 
only covered a 83 cm long section, on the central line of the model, 
and allowed us to image up to 18 cm of depth. 


4.2.2. Construction of BBAs from geophysical and geotechnical data 
The electrical acquisition was carried out on 83 cm long, on the cen- 
tral line of the model, with a first acquisition on 47 cm, and three next 








Fig. 14. Transparency view of the test bench. 


Table 1 
Values of electrical resistivity and depth of penetration of the materials set up within the 
test bench. 

Plaster before pluviation Saturated Hostun sands 


Electrical resistivity (Q.m) 


Mean 31.28 78.15 
Standard deviation 3.23 11.18 
Number of measures 12 52 


Penetration depth (mm) 


Mean 0.11 17.31 
Standard deviation 0.04 1.61 
Number of measurements 8 10 
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Fig. 15. Scheme of the idealized section model (with vertical exaggeration), including the FoD constituent events associated with the materials of the test bench. 


acquisitions done after respective displacements of 12 cm (roll along 
method). The results obtained from the inversion of the acquired data 
are displayed in Fig. 16a. These results make it possible to highlight 
the existence of three distinct sets, at depths relatively close to the tar- 
get model (Fig. 15) but presenting vertically slightly shifted interfaces, 
gradual rather than sharp. In addition, the variation in saturated sand 
thickness is poorly evaluated. Indeed, the anomalous zone is recognized 
but associated here, in its lower part, with values of electrical resistivi- 
ties much larger than what they really are. 

The proposed values, although in the same order of magnitude, do 
not exactly match the ranges of values measured on the materials inde- 
pendently (Table 1). In order to characterize the events (materials) of 
the FoD, a distribution in modal classes (Fig. 16b) is used to visualize 
the number of cells of the discretized section for the 2D inversion, asso- 
ciated with their corresponding ranges of resistivities. This distribution 
makes it possible to highlight the three large sets of materials in the 
model. Thanks to it, the bounds of the events considered can thus be de- 
fined, in Q.m, so that the intervals are the same length, as presented 
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Eq. (18): 


(10; 35] 

(40; 140] 

[9 500; 33 250] 

(2.85; 10[U]35; 40[U]140; 9 500[U]33 250; 116 375] 


1 


(18) 


0 
02 
03 
04 


In contrast to information from the geophysical source, geotechnical 
data were obtained beforehand by laboratory penetration cone testing, 
and then numerically simulated prior to fusion. Several simulations pro- 
posing various positions of survey points were carried out. In order to 
simulate drilling points, the associated mean depth values (mm) and as- 
sociated standard deviations (Table 1) were used to draw values, fol- 
lowing a normal distribution defined for each event. An average 
penetration depth value of 0 mm is used for 03 (bottom of tank) and 
an associated standard deviation of 0.01 mm, meaning that negative 
values may be drawn. These random draws, make it possible to define 
the limits, in mm, of the intervals associated with the events of the 
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Fig. 16. a) inverse model resistivity section obtained by roll along acquisitions in the central line of the model and b) model classes’ distribution of the cells presented in a), according to the 


electrical resistivity values (Q.m). 


12 T. Dezert et al. / Journal of Applied Geophysics 170 (2019) 103824 


FoD as presented Eq. (19): 


6, = 0.04; 0.19] 

6, = [13; 21] 

6; = [—0.02; 0.02] 

64 = [—0.05; —0.02[U]0.02; 0.04[U}0.19; 13[U]21; 100} 


(19) 





Thus, 2 mm wide boreholes are simulated, down to 15 cm and ac- 
quiring every 5 mm with an associated inaccuracy of 0.01 mm. The 
values of penetration depth obtained can then be associated with the 
different materials of the model. 


5. Test bench data fusion results 


The results of the geophysical and geotechnical information fusion, 
are proposed in Fig. 17. The simulations were carried out according to 
four distinct vertical drill positioning configurations, represented in 
dashed lines in the figures and at regular intervals: i) 8 holes inter- 
spaced of 10 cm (Fig. 17i) (x = 10; 20; 30; 40; 50; 60; 70; 80 cm), ii) 
5 holes inter-spaced of 18 cm (Fig. 171i) (x = 4, 22, 40, 58, 76 cm), iii) 
3 holes inter-spaced of 25 cm (Fig. 17ili) (x = 15, 40, 65 cm), iv) 2 
holes inter-spaced of 50 cm (Fig. 17iv) (x = 15, 65 cm). The fusion re- 
sults carried out are presented, respecting i) the hypothesis of Smets 
(Fig. 17a and b, ii) the hypothesis of a closed-world (Section 2.1) with 


Considered events 





PCR6 rule (Fig. 17c and d). Fig. 17b and d represent the belief mass 
values associated with events having the largest mass, represented re- 
spectively in Fig. 17a and c. The fusion results are analyzed and 
discussed in the next section. 


6. Fusion results analysis and discussion 
6.1. Different rules of combinations 


Let us discuss and compare the results obtained by the 2 different 
combination rules used in an 8-boreholes simulation (Fig. 17i). In the 
framework of a model as rich in geotechnical information, the section 
proposed by the PCR6 method (Fig. 17i.c) is very close to the target 
model set up (Fig. 15). The three sets are well characterized and the in- 
terfaces at 2.5 cm deep (sands-plaster) and at 7.5 cm deep (plaster-PVC 
tank and sand-PVC tank) are much better defined than by ERT alone 
(Fig. 16a). Moreover, thanks to this geotechnical information, the sand 
thickness anomaly could be correctly characterized as saturated sands 
(02) and not as a more resistive anomaly, in continuity with the insulat- 
ing material from below, as suggested by the results of the inversion. 
The lateral extension of this anomaly is, moreover, well estimated 
(20 cm). The combination of Smets highlights the significant conflict 
existing between the two considered sources of information (Fig. 171. 
a) concerning the two first layers. 
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Fig. 17. Representation of the events having the highest belief mass in a) and c) and their associated mass values in b) and d). For i) 8 boreholes, ii) 5 boreholes, iii) 3 boreholes, iv) 2 
boreholes are considered. For each case, (a,b) figures are results of Smets fusion, while (c,d) figures are results of PCR6 fusion. The borehole positions are in dashed lines. 
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Whatever method is used, the presence of a hypothesis 0, is found at 
the vertical and horizontal interfaces (Fig. 17i.a, c). This hypothesis does 
not correspond to any material set up in the test bench. The belief 
masses attributed to such a hypothesis, highlight the transition zones 
not conform to reality, proposed by the inversion of the electrical resis- 
tivity data (Fig. 16a). In comparison to the belief masses associated with 
the other hypotheses of the model, the belief masses associated with 04 
are the lowest (Fig. 171.b, d), showing that the confidence granted to 
such a material remains quite relative. An overall confidence drop is 
also observed from 15 cm depth. This corresponds to the maximum 
depth reached by the simulated boreholes. As confidence is extended 
laterally, the belief masses are constrained only by geophysical informa- 
tion to such a depth and therefore rely only on one source of 
information. 


6.2. Influence of the number of boreholes and positions 


The first intuition would be to assume that the more the number of 
boreholes decreases, the more the method should be put in difficulty to 
properly characterize the section of the set up test bench. Although this 
is partly true, the quality of the results is not based as much on the num- 
ber as on the positions of the drillings. Indeed, the anomaly of saturated 
sands contained between the two banks of plaster (Fig. 15) is as well 
characterized in terms of lateral extension with three or five soundings 
(Fig. 1 7ii.c, ili.c). It also has an equivalent associated trust (Fig. 17ii.d, iii. 
d). It turns out that the belief masses associated with the event 0; (plas- 
ter) are even smaller for a fusion including three soundings (Fig. 1 7iii.d) 
than for a simulation of only two (Fig. 17iv.d). 

The explanation of such results lies in the fact that being in the pres- 
ence of consecutive boreholes, informing about the occurrence of differ- 
ent materials, at an equivalent depth, induces a rapid decrease in the 
confidence attributed to the boreholes. Therefore, more credibility is 
given to the geophysical information source, explaining the greater 
presence of 04, which reflects the gradual transitions in electrical resis- 
tivities. The masses associated with this event, however, remain rela- 
tively small. On the other hand, if two consecutive boreholes have the 
Same geotechnical values, for a specific depth, the lateral decay rate 
will be low and no priority can be given to a different material existing 
between these two boreholes. That is why the sand anomaly in the cen- 
ter of the model does not appear in the results fusion with two sound- 
ings (Fig. 17iv.c): no borehole pass through the anomaly and the 
geophysical source is unable to characterize this material as saturated 
sand. The strength of these results is that they suggest the presence of 
04 in this location, suggesting that the survey campaign should be rein- 
forced (with a new borehole position for example). 

The conflict presented by Smets combination (results in Fig. 17i.a, ii. 
a, lil.a and iv.a) is neither a function of the number of geotechnical 
soundings. In this study, the cases of fusion bringing the highest amount 
of conflict are in fact the ones with eight and two soundings. Nor is it to 
be confused with a lack of knowledge of the subsoil. Conflict zones high- 
light contradictory information between the two sources. These zones 
are generally between two consecutive boreholes providing the same 
information, but going against the available geophysical information. 
These are therefore potentially anomalous zones where the geophysical 
information must be considered carefully, in particular if the belief mass 
associated with the event retained after normalization is too low. 


6.3. Important considerations and potential in the application 


It is important to consider that the effectiveness of this fusion meth- 
odology has been assessed by comparing the fusion results with a target 
model (Fig. 15). However, this remains an idealized representation of 
the test bench set up and could be, in some places, quite far from reality 
(real interfaces not perfectly horizontal or vertical, materials not per- 
fectly homogeneous, 3D effects neglected ...). The approach is different 
from the one of the synthetic study (Fig. 4) where the model shown 


corresponds to the true model. In order to control the effectiveness of 
the fusion methodology, it was envisaged to carry out ex post verifica- 
tions of the constituent materials. Unfortunately, for practical reasons, 
this could not be done (reworking of materials modifying their physical 
properties, interaction with water, delicate cutting and extraction ...). 

Regarding the fusion methodology developed, two aspects are de- 
batable. First, the choice to set a mass of belief equal to 1 on the geotech- 
nical information in the boreholes. Second, the effect of different 
random draw results on the fusion results. The choice of a maximum 
punctual confidence (m = 1) in boreholes is defended in order to give 
a full and local confidence to geotechnical information as it is currently 
done during investigation campaigns. Excessive risks are not taken since 
the test bench is relatively well known and the synthetic model is per- 
fectly well known. Thus, it is sure that simulated borehole values refer 
to the right materials. Furthermore, a value of m = 0.99 instead of 
m = 1, for instance, does not significantly change the results and does 
not change the interpretation and the resulting discussion. Regarding 
the effect of random draws, these draws were done following a normal 
distribution, the variations from one draw to another are minimal and 
the results of fusion differ little. 

Such an information fusion algorithm, dedicated to the combination 
of data from geophysical and geotechnical sources, should prove useful 
for processing of data acquired during investigation campaigns for 
many different kinds of issues. It is possible to envisage its use with a 
larger number of materials, but also, and especially, with a larger num- 
ber of data types from geophysical methods (seismics, ground penetrat- 
ing radar) and geotechnical testing methods (penetration cone, core 
Sampling with laboratory identification, permeability tests ...) 
associated. 

In the framework of a recognition campaign, the conflict zones, or 
zones with a low associated confidence, would make it possible to spec- 
ify the locations where the investigation must be reinforced. The ulti- 
mate goal is to obtain a more robust and cost-effective diagnosis of 
the investigated structure, more targeted for geotechnical investigation. 
This methodology has particularly shown its ability to correctly charac- 
terize interfaces, which corresponds to areas where the risk of instabil- 
ity is potentially the greatest. For a levee embankment issue, for 
example, the results from such a methodology could come to feed into 
models of breakage risks (ex: CARDigues (Apel et al., 2004) ). 


7. Conclusion 


In this work, a new methodology has been presented, based on belief 
functions to take benefit and to combine two different and complemen- 
tary kinds of information: geophysical and geotechnical. Each one 
having its own spatial distribution and related uncertainties and inaccu- 
racies. A new representation of the information has been proposed, tak- 
ing into consideration two different investigation methods, associated 
with degrees of belief. This representation is more informative than 
data superposition of different physical parameters. 

In the first place, this new approach has been validated with a syn- 
thetic study, simulating data acquired by ERT and a CPT method, consid- 
ering a 2D model with two layers and thickness variation. The results 
were obtained with different noise ratios applied to the geophysical 
data and different values of lateral decay coefficient for the geotechnical 
information. The most appropriate value to pick up for the coefficient 
has been pointed out and it has been showed that this approach 
was able to manage the noise ratio, thanks to the use of Wasserstein 
distances. 

In order to address the problem of combining information acquired 
by geophysical and geotechnical methods during investigation cam- 
paigns, and to acquire values from real devices, a test bench composed 
of plaster and saturated sands was set up. The methods used to charac- 
terize such a physical model were the ERT method (geophysical) and 
the laboratory penetration cone method (geotechnical). While the 
data has been acquired by a dedicated small scale ERT device, on the 
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surface and on the central line of the complete model, borehole 
were simulated respecting the penetration depth ranges previously 
established. 

Fusion results were proposed following 2 combination rules (Smets 
and PCR6) as well as for four different simulations of number and posi- 
tions of boreholes. The results highlighted the ability of this fusion ap- 
proach to correctly characterize the test bench materials as well as to 
specify the positions of the interfaces (vertical and horizontal) between 
the materials. Moreover, for each result, thanks to a graphical represen- 
tation, the associated confidence is proposed. 

Further research should include cases of material mixtures and cases 
of different materials sharing common ranges of physical properties in 
order to test the ability of this methodology to differentiate them. We 
also wish to test this new methodology in real investigation campaigns 
in order to improve the available knowledge and strengthen the charac- 
terization. The level of confidence associated with the proposed results 
may be very relevant for decision support (eg models of failure haz- 
ards). The results of such a methodology should make it possible to pro- 
pose the most relevant borehole positions (that are a function of 
conflictual and anomalous areas), in order to make the quality of the in- 
formation more cost-effective. 
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